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Abstract. We evaluate the equilibrium state of a mixture of "^Li and ^Li atoms with 
repulsive interactions, confined inside a pancake-shaped trap under conditions such 
that the thickness of the bosonic and fermionic clouds is approaching the values of the 
s-wave scattering lengths. In this regime the effective couplings depend on the axial 
confinement and full demixing can become observable by merely squeezing the trap, 
without enhancing the scattering lengths through recourse to a Feshbach resonance. 
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1. Introduction 

Dilute boson-fermion mixtures have recently been studied in several experiments by 
trapping and cooling gases of mixed alkali-atom isotopes [H 12 13 IH 13 El • The boson- 
fermion coupling strongly affects the equilibrium properties of the mixture and can lead 
to quantum phase transitions. In particular, it has been shown that boson-fermion 
repulsions in three-dimensional (3D) clouds can induce spatial demixing of the two 
components when the interaction energy overcomes the kinetic and confinement energies 
[7j. Several configurations with different topology are possible for a demixed cloud 
inside a harmonic trap Although spatial demixing has not yet been experimentally 
observed, the experiments of Schreck et al. P on a ^Li-^Li mixtures inside an elongated 
3D trap appear to be not far from the onset of a demixed state: an increase in the boson- 
fermion scattering length by a factor five would be needed to enter the phase-separated 
regime [H]. 

In this Letter we examine the possibility of attaining phase separation by merely 
varying the trap geometry for the ^Li-''Li mixture at the "natural" values of the 
scattering lengths. We focus on the case of quasi-two-dimensional (Q2D) confinement 
inside a pancake-shaped trap. We find that this geometry favours the demixed state 
through the appearance of the axial size of the atomic clouds in their effective coupling 
in the azimuthal plane. We also show that several configurations are possible in the 
plane of the trap at zero temperature, but the configuration consisting of a core of 
bosons surrounded by a ring of fermions is the most energetically favourable for a wide 
range of values of system parameters. 



2. The model 

The atomic clouds are trapped in pancake-shaped potentials given by 

^7 (x, z) = mbjujl^,j){x^ + \ljy^)/2 + mbjuj\^j~^z^ /2 = Vbj{x, y) + Ubj{z) (1) 

where m^j are the atomic masses and ujz{b,f) ^ ^x{bj) the trap frequencies. We 
focus on the case where the trap is fiat enough that the dimensions of the clouds 
in the axial direction, which are of the order of the axial harmonic-oscillator lengths 
(^z{bj) = {fi'/^bj^z{bj)Y^'^, comparable to the 3D boson-boson and boson-fermion 
scattering lengths abb and a^f. In this regime we can study the equilibrium properties 
of the mixture at essentially zero temperature (T ~ 0.02 Tp) in terms of the particle 
density profiles in the {x,y} plane, which are nc{x,y) for the Bose-Einstein condensate 
and nf{x,y) for the fermions. The profiles are evaluated using the Thomas- Fermi 
approximation for the condensate and the Hartree-Fock approximation for the fermions. 
We shall take a^t = a^f (= a^, say) and only near the end we shall discuss the case 

The Thomas-Fermi approximation assumes that the number of condensed bosons 
is large enough that the kinetic energy term in the Gross-Pitaevskii equation may be 
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neglected [TU]. It yields 

nd^^y) = [fJ'b - Vb{x,y) - gbfnf{x,y)]/gbb (2) 

for positive values of the function in the square brackets and zero otherwise. Here, 
/if, is the chemical potential of the bosons. This mean-field model is valid when the 
high-diluteness condition ncfl^^ -C 1 holds and the temperature is outside the critical 
region. If the conditions > abb,abf are fulfilled, the atoms experience collisions in 
three dimensions and the coupling constants in Eq. Q can be written in terms of the 
3D ones as jlT] 

„3D „3D 

^ _ 9bb ^ _ 9bf 

9bb - r— , 9bf - r— \P) 



with gf^ = 4:71 abb /iTib, 9bf — 'i.nt^abf /rrir and rrir = mbmf/{mb + T^f)- 

The Hartree-Fock approximation [T^ IT^ treats the fermion cloud as an ideal 
gas subject to an effective external potential, that is 



nf{x,y) 



(4) 



where /x/ is the chemical potential of the fermions and 

^^//(x, y) = Vf{x, y) + gbfUc^x, y). (5) 

The fermionic component has been taken as a dilute spin-polarized gas, for which the 
fermion-fermion s-wave scattering processes are inhibited by the Pauli principle and 
p-wave scattering is negligibly small ^Sj- In the mixed regime the Fermi wave number 
in the azimuthal plane should be smaller than l/of,/, but this is not a constraint in the 
regime of demixing where the boson-fermion overlap is rapidly dropping. 

The chemical potentials iXbj characterize the system in the grand-canonical 
ensemble and are determined by requiring that the integrals of the densities over 
the {x, y} plane should be equal to the average numbers Nb and Nf of particles. 
The presence of a bosonic thermal cloud can be treated by a similar Hartree-Fock 
approximation P|, but it has quite negligible effects at the temperatures of present 
interest. 



3. General conditions for phase separation 

The effective strength of the atom-atom collisions in the azimuthal plane depends in 
our Q2D model on the axial harmonic-oscillator length according to Eq. 0. We derive 
and illustrate in this section the consequences for demixing at very low temperatures. 

We preliminarly recall that in the macroscopic limit the phase transition is sharp 
and the overlap between the two species after demixing is restricted to the interfacial 
region, where it is governed by the surface kinetic energy (see Ref. jTH] for an example 
in a 3D model). In the case of mesoscopic clouds under confinement the transition is 
instead spread out and several alternative definitions of the location of demixing can 
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therefore be given. We discuss below three alternative locations, that we denote as 
partial, dynamical, and full demixing. Their definition and the derivation of simple 
analytical expressions in the Q2D model are given in the following. 

3.1. Partial demixing 

The interaction energy Eint between the boson and fermion clouds initially grows as the 
boson-fermion coupling is increased, but reaches a maximum and then falls off as the 
overlap between the two clouds diminishes. We locate partial demixing at the maximum 
of Eint, that we calculate from the density profiles according to the expression 



In the calculations that we report in this section we have used values of system 
parameters as appropriate to the experiments in Paris on the ^Li-^Li mixture P, that 
is Nk = NfC^ W and cj^fe/27r = 4000 Hz, u;^//27r = 3520 Hz. 

We first fix the two scattering lengths at their "natural" values (a;,;, = 5.1 Oq and 
flfc/ = 38 Oo, with oo the Bohr radius) and vary the thickness of the trap from 100 a;,/ 
down to abf. By this choice we are increasing both the boson-boson and the boson- 
fermion effective coupling in the azimuthal plane. The behaviour of the interaction 
energy as a function of the ratio a^j jaz at fixed abb and abf is shown in Fig. [H It is seen 
that Eint goes through a maximum at ~ 16 Ob/. 




(6) 
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Figure 1. Boson-fermion interaction energy Ei„t in a circular disc (in units of Nhuj^b, 
witli N = Ni) + Nf), as a function of ai,f /a^ for fixed values of ati, and a^f. The arrows 
indicate the locations of demixing estimate from Eqs. |7J, (|12|l and (|14(l . 

As an alternative and with the aim of a later comparison with the results previously 
obtained in the 3D case [HI El 5 we show in Fig. |21the interaction energy as a function of 
the ratio a(,//abf, at fixed a^;,. The axial thickness has been set equal to the largest of 
the two scattering lengths. This condition fixes the limit of validity of our model jl7j . 
The maximum of Eint lies at abj/ahi, ~ 0.8, much below the value abf/abb ~ 7 of the 
ratio of the "natural" scattering lengths. In fact, at abf/abb — 7 the two clouds would 
undergo demixing as the axial trap stifness is increased towards a Q2D model. 
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Figure 2. Boson-fermion interaction energy Eint in a circular disc (in units of Nhujxb), 
as a function of a^f /abb for abb = 5.1 ag. The arrows indicate the locations of demixing 
estimate from Eqs. 0, and (^Jj). 



An approximate analytic formula for the location of partial demixing as a function 
of the system parameters in the Q2D model can be obtained from the condition 
dEint/dghf = by using the estimate Ucj ~ N^j / {T.tiRIj) with the values of the 
cloud radii in the absence of boson-fermion interactions, Rj = (SNj / XfY^^a^f and 
Rb = {16Nhabb/ {^/27Taz\b)y^'^axb where ax(bj) = (fi/'mbj^xibj))^^'^- The location of the 
maximum in the boson-fermion interaction energy as a function of atf/abb lies at 

1/2 



where 



with 



and 



CLbb 



part 



Ipart ( — ) (7) 



a-bb 




Ipart = Ci-j-j^ + 02— rj^ (8) 



c.={'-Y^^(^Y^ (9) 



TT / 2mrUJxb V 



c2={^Y'^^^(^y'\ (10) 
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We recognise in Eq. (jHl) a geometric combination of two scaling parameters: ciNj / nJ 

1 /2 1 /2 

is dominant in the case A*";, <^ Nf while C2N^ /Nj is dominant in the opposit limit. 
The prediction obtained from Eq. (jZj) is shown in Figs. Hand El by a solid arrow. There 
clearly is good agreement between the analytical estimate and the numerical results. 



3.2. Dynamical demixing 

On further increasing the boson-fermion coupling one reaches the point where the 
fermion density vanishes at the centre of the trap. This occurs when 

^ (11) 

9bb dyn l^b 

We denote this point as the dynamical location of the demixing in a mesoscopic cloud, 
since we expect a sharp upturn of the low-lying fermion-like collective mode frequencies 
to occur here as it is the case for both coUisional and coUisionless excitations in a mixture 
under 3D confinement IT^ . 

If we insert the chemical potentials for ideal-gas clouds in Eq. flllll . we obtain an 
approximate expression for the location of dynamical demixing, 

1/2 



dbb 

where 



dyn 



Idyn [ ^ ] (12) 
O'bb 



,,^^nrf^^^{ML)"\ (13) 

\2J mb + nif \NbXb J 
The condition for dynamical demixing coincides with that for partial demixing in the 
limit Nb ^ Nf. The prediction obtained from Eq. ()12|) is indicated in Figs. ^ and |21 by 
a dot-dashed arrow. 



3.3. Full demixing 

The point of full demixing is reached when the boson-fermion overlap becomes negligible 
as in a macroscopic cloud. Using the instability criterion outlined in Ref. jSHI for the 
3D case, the condition for full phase separation at T = is 

1/2 



abf_ 

Clbb 



full 



where 



Vuiii^] (14) 

7/., = / . (15) 

At variance for the condition for full phase separation in 3D, for the Q2D confinement 
this criterion does not depend on the number of fermions, while l/a^ plays the role of 
the Fermi momentum. 
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The prediction obtained from Eq. ((Tj) is indicated in Figs. [T]and|2]by a short-dashed 
arrow. 

3.4- Density profiles at full demixing 

In summary, it can be seen from Eqs. ((Zj), (|T2|l and ()14|1 that the two control parameters 
for the transition in a mesoscopic cloud under Q2D confinement are abf/obb and az/abb- 
Our main conclusion is that in a Q2D geometry a relatively small value of the boson- 
fermion scattering length suffices to reach the regime of full phase separation. In the case 
of the ^Li-^Li mixture the bare values of the boson-boson and boson-fermion scattering 
lengths fulfill the condition for full phase separation for values of lower than 10 a^/. 

In view of the above result, it is useful to illustrate in Fig. 01 the density profiles 
to be expected in the fully demixed regime for a mixture under isotropic confinement 
in the azimuthal plane (notice that a circular disc appears in Fig. IHlas an oval, because 
of the different scales used on the horizontal and vertical axes). In fact we have found 
various metastable configurations for the demixed mesoscopic cloud in addition to the 
thermodynamically stable one consisting of a core of bosons surrounded by a ring of 
fermions. 

In Fig. 3 we show topviews of the density profiles in the {x, y} plane for the choice 
of parameters 

(^bb — 5.1 flg, a,bf — 38 flg and = Obj. The most energetically stable 
configuration is shown in Fig. Efa) and lies at energy E = 74.90 Nhwxb- Figure Efb) 
shows a configuration with fermions at the centre surrounded by a ring of bosons inside 
a fermion cloud {E = 78.94 Nhwrcb)- In Fig- Efc) a fermion slice lies between two boson 
slices surrounded by a fermion cloud {E = 79.04 NfkJj.b). Finally, Fig. ^d) shows an 
asymmetric configuration in which a core of bosons is displaced away from the centre of 
the trap {E = 80.20 Nhw^b). The latter two configurations break the symmetry of the 
trap and this is an effect due to the finite size of the system. 

Configurations (b)-(d) are structurally but not energetically stable, i.e. they 
represent metastable structures in local energy minima. They have been obtained by 
using density profiles with various topologies as the initial guess for the self-consistent 
solution of Eqs. Q and (jH). In experiments bulk modes with the appropriate symmetry 
may be exploited to attain these "exotic" configurations pij . 

We have also studied the behaviour of the configurations shown in Fig. IHl as a 
function of the scattering lengths (S.Ioq < abb < S x IO^Cq and 38ao < a?,/ < 5 x 10 ''oo) 
and the planar-anisotropy parameters (10~^ < Xbj < 400). In contrast to the case of 
an elongated 3D confinement studied in our previous work [S], where the lowest-energy 
configuration can have different topologies depending on the system parameters, we 
have found that the configuration shown in Fig. ISfa) remains the most energetically 
stable for the Q2D geometry. A configuration consisting of a central core of fermions 
surrounded by a cloud of bosons may, however, be the stable one for sufficiently strong 
boson-boson repulsions [S]. 
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Figure 3. Density profiles at full demixing for fermions (left) and bosons (right) in 
the {x,y} plane, for a Q2D mixture with abb — 5.1 oq, Ob/ = 38 ao, and \b — \f = 1: 
(a) core of bosons surrounded by a ring of fermions; (b) core of fermions surrounded 
by a ring of bosons inside a fermion cloud; (c) fcrmion slice between two boson slices 
surrounded by a fermion cloud; (d) core of bosons displaced away from the centre of 
the trap and surrounded by a fermion cloud. Notice the difference in horizontal and 
vertical scales. 
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4. Summary and concluding remarks 

In summary, we have focussed on a boson-fermion mixture confined inside pancake- 
shaped traps, such that the scattering events can still be considered three-dimensional 
but nevertheless affected by the vertical confinement. By using a mean- field description 
for the equilibrium densities in the azimuthal plane, we have studied the boson-fermion 
interaction energy as a function of the thickness of the atomic clouds and of the boson- 
boson and boson-fermion scattering lengths. We have given approximate analytical 
expressions identifying three critical regimes in terms of the scaling parameters afe//afeb 
and az/abb'- (i) partial demixing where the boson-fermion interaction energy attains 
maximum value from a balance between increasing interactions and diminishing overlap; 
(ii) dynamical demixing where the fermionic density drops to zero at the centre of 
the trap and a sharp dynamical signature of demixing may be expected; and (iii) full 
demixing where the boson-fermion overlap is negligible as in the macroscopic limit. 
These different criteria for the location of the phase transition in a mesoscopic cloud 
yield quite close predictions in a pancake-shaped trap under strong axial confinement. 

We have found several metastable configurations for the demixed cloud, having 
various topologies but lying at higher energy above the stable configuration which is 
composed by a core of condensed bosons surrounded by fermions. We have verified 
that this is the minimum-energy configuration over extensive ranges of values for the 
coupling constants, the anisotropy in the plane of the trap, and the thickness of the 
trap. 

A main result of our work is that full demixing can be reached in the ^Li-^Li 
mixture by merely tuning the thickness of the trap, without necessarily tuning the 
scattering lengths. The full-demixing critical value of about 10 a?,/ for the thickness is not 
attainable in the actual experiments with the bare value of the ^Li-^Li scattering length, 
and a combination of squeezing of the pancake-shaped trap and exploiting Feshbach 
resonances should be used. On the contrary, the dynamical location of demixing may 
be observed by only varying the number of particles, without enhancing the scattering 
length. In fact for a mixture of 10^ atoms of ^Li and 10^ atoms of ^Li trapped in the 
radial confinement discussed in Sec. 13.11 it may be possible to observe a sharp upturn of 
the low-lying fermionic modes for a pancake thickness of the order of ~ 600 abf — 1 yum. 
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acknowledges support from TUBITAK and from the Research Fund of the University 
of Istanbul under Project Number 161/15102003. 
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